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Abstract 

In this paper we study the spectrum of the spin-1 Temperley-Lieb 
spin chain with integrable open boundary conditions. We obtain the 
eigenvalue expressions as well as its associated Bethe ansatz equations 
by means of the coordinate Bethe ansatz. These equations provide the 
complete description of the spectrum of the model. 
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1 Introduction 



Quantum integrable models and their associated classical vertex models have 
been largely investigated over the years [TJ [2] . Many of these models have 
been widely studied by Bethe ansatz techniques with different boundary con- 
ditions and at zero or finite temperature and magnetic field. 

Nevertheless one still has interesting problems withstanding to the stan- 
dard techniques. One of these problems is the biquadratic spin-1 model 
[21HJ|5]. The biquadratic model was shown to be invariant by the Temperley- 
Lieb algebra |6J . This property has made possible the discovery of new inte- 
grable quantum spin chains, which was achieved by exploiting the represen- 
tation theory of the Temperley-Lieb algebra [TJ [8] . 

Moreover the biquadratic model and its generalizations (Temperley-Lieb 
spin chains) were solved by coordinate Bethe ansatz for periodic boundary 
conditions and free ends [31 El [101 E]. Afterwards, it has also appeared 
the spectrum of the transfer matrices by means of functional methods [12]. 
However, there is still no algebraic Bethe ansatz formulation for these models. 

Recently the concept of Temperley-Lieb equivalence [T] was used in or- 
der to obtain the spectral properties of quantum spin chains of Temperley- 
Lieb type for periodic boundary conditions and free ends [13]. The study of 
the spectral multiplicities allowed to the computation of the thermodynamic 
properties at finite temperature [T3] . 

Besides that, the solution of reflection equation, associated to the problem 
of Temperley-Lieb spin chains with integrable open boundaries, was recently 
obtained [HI [T5]. However, the computation of the spectra of these spin 
chains is still an open problem. 
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In order to fill this gap, in this paper we are interested in the spectra 
of the U q [sl(2)] Temperley-Lieb spin chain with integrable open boundary 
conditions [H]. We use a suitable generalization of the coordinate Bethe 
ansatzfTU] in order to obtain the eigenvalues of the spin-1 Temperley-Lieb 
model with diagonal open boundaries. 

The outline of the article is as follows. In section |2] we introduce the 
Temperley-Lieb spin chain with integrable open boundaries. In section [3] 
we discuss the application of the coordinate Bethe ansatz and obtain the 
eigenvalues of the Temperley-Lieb spin chain. Our conclusions are given in 
the section HI 

2 Temperley-Lieb spin chain 

The cornerstone of the theory of quantum integrable models in one-dimension 
is given by the Yang-Baxter equation, 



This equation provides the commutativity property of the transfer matrix 



T(A) = Tr.4 [T A (\)}, where T A (X) = R AL (X) ■ ■ ■ R A1 (X) , R 12 (X) = P 12 R 12 {X) 



and P\2 is the permutation operator. Quantum integrable spin chain with pe- 
riodic boundary conditions are obtained by means of the logarithmic deriva- 
tive of the transfer matrix T(A). 

The Temperley-Lieb invariant solutions of the Yang-Baxter equation 
are well known[16j. The spin-1 U q [si (2)] solution can be written as, 



Ru(X - ^R 23 (X)RM = R 2 MR 12 (X)R 23 (X - //). 



(1) 



Rij(X) 



sinh(7 — A) 



sinh A 



Uij, 



(2) 



sinh 7 



sinh 7 
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where 2 cosh 7 = g 2 + 1 + g 2 and the Temperley-Lieb operator is given by 
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(3) 



which is the projector onto the two-sites spin zero singlet written in the basis 
{|+> , |0> , |-». 

The notion of integrability was extended to tackle with open boundary 
problems [T7]. On the one hand, the .R-matrix describes the bulk dynamics, 
and on the other hand a new set of matrices, the if-matrices, represent the 
interactions at the left and right ends of the open spin chain. This is a 
consequence of the reflection equation, which reads 

R 12 {\- H)K 1 (\)R 21 (\ + n)K 2 (ti) =K 2 (ji)R 12 {X + f i)K 1 {X)R 21 {X-n). (4) 



In the case of open boundary conditions, the transfer matrix can be writ- 
ten as 

f(A) = Tt A [^ +) (A)T^(A)^- ) (A) [^(-A)]- 1 ] , (5) 

where K^(X) can be chosen as one of the solutions of the reflection equation 
(4). The other boundary matrix K^(X) is obtained from the previous one 
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by means of the isomorphism 

K { 2\X)=K ( X \-\-p) t V% 



(6) 



where t means transposition, the crossing parameter is p = —7 and the 
crossing matrix is given by 

/ n n „\ 



V 









Q 





-1 





-1 









(7) 



/ 



The integrable open spin chain is obtained by means of the logarithmic 
derivative of the transfer matrix (pi), such that, 

simi7 d 



H 



2 dX 

L-l 

/] Uk,k+i + 



lnt(A) I A + const, 
sinhj dK[~\\) 



(8) 



k=l 



d\ 



A=0 



+ 



Tr^ 


~k£\o)u la 


Tr 


A 







The solutions of the reflection equation (|4]) associated to the i?-matrix 
(|2]) were recently obtained [H]. Here we list only the diagonal ones, which 
we will use throughout this work. 



^(1,0,0) (*) 



( k [ (l(\) \ 
1 
1 



^(0,1,0) (*) 



^(0,0,1) ( A ) 



( 1 N 

kg (A) 
II 



1 
1 
V 4 7 3 ] (A) J 



(9) 
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where 



*S(A) 

4!i(A) 
and 



#c 2 (A) [(1 + g 2 )x 2 (A) + xx(A)] + 2 MA)x' 2 (A) - x > 1 (X)x 2 {\)] 
-Px 2 (\) [q~ 2 x 2 (\) + x 1 (X)} + 2 MAy 2 (A) - xi(A)x 2 (A)] ' 
(3x 2 {\) [(g- 2 + g 2 )x 2 (A) + ^(A)] + 2 [x 1 (A)x 2 (A) - x' 1 (X)x 2 (\)] 

-(5x 2 (X) [x 2 (X) + x l (X)} + 2 [x 1 (A)x 2 (A) - x' 1 (A)x 2 (A)] 
(3x 2 {X) [(q- 2 + l)x 2 (X) + xx(X)\ + 2 [xi(A)4(A) - x> 



-f3x 2 (X) [q 2 x 2 (X) + Xl (X)\ + 2 MA)x' 2 (A) - x[(X)x 2 (X)\ 
sinh(7 — A) , , s sinh A 



an(A) 



.^(A) 



(11) 



sinh 7 sinh 7 

There is an additional set of diagonal solutions of different form, given by 



K 



(-) 
(0,1,1) 



(A) 



( X 

fcf?(A) 











(-) 

(1,1,0) 



(A) 



^3,3 W 





\ 



*£?(A) y 

N 



(-) 

(1,0,1) 



(A) 



4,2 (A) 













k 



[ii] 
2,2 



(A) ) 



45(A) 







(12) 



where 



45(A) 



45(A) 



45(A) 



f3x 2 (X) [g- 2 x 2 {X) + xi(A)] + 2 MA)x 2 (A) - x / 1 (A)x 2 (A)] 
-f3x 2 {X) [(1 + g 2 )x 2 (A) + x 1 {X)} + 2 [a;i(A)4(A) - x;(A)x 2 (A)] ' 

Px 2 {X) [x 2 (X) + Xl (X)] + 2 [xx(A)x 2 (A) - x;(A)x 2 (A)] 
-/3x 2 (X) [(q- 2 + g 2 )x 2 (A) + ari(A)] + 2 [m(AK 2 (A) - x / 1 (A)a: 2 (A)] 
/3x 2 (A) [g 2 x 2 (A) + 2* (A)] + 2 [xi(A)x 2 (A) - x , 1 (X)x 2 (X)] 



-f3x 2 (X) [(q- 2 + l)x 2 (A) + Xi(A)] + 2 [xi(A)x / 2 (A) - x[(X)x 2 

The boundary terms of Hamiltonian ^ are directly obtained from the 
above K^' matrices. In particular, the left boundary acting non-trivially in 
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the site 1 has the form 



Bi 



sinh 7 dK { ~\\) 



dX 



A=0 



V 



hi o o 

l 22 
l 33 



while the right boundary acting in the last site L is given by 



B, 



Tr^ 


K ( +\0)U L>A 




~k a + \o) 





( rn ^ 

r 22 

\ o o r 33 y 



where 



^33 



9" 2 *£(0) 



fc+ (0) + *+(0) + A;+(0) 



r 2 2 



*&(0) 



fc+(0) + /c+(0) + /c 3 + 3 (0)' 



k+(0) + fc+(0) + fc 3 + 3 (0)' 
Therefore, we can compute the left boundary terms (14) as 

/ 



B\ 



(3,m) 



V 



/3 sinh 7 
2 
































P sinh 7 
2 



\ 



/l 



5} 



(2,m) 



ooo 






/l 



/l 



corresponding to the left i^-matrices i^l ^ Q w f^ l -[ ^ and ^ ^ respectively. 
Similarly, we have three more left boundaries corresponding to the K- 



(14) 



(15) 



(16) 



(17) 
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matrices ^(0,^1)^(1,0,1) and #(1,1,0) 



(urn 



\ 





g /3 sinh 7 





2 





' /3 sinh 7 
2 



(3,[//]) 







/3 sinh 7 
2 





(2,[JJ]) 



^ /3 sinh 7 
2 















/9 sinh 7 



V 








/3 sinh 7 q 
2 U 





where /3 is the left boundary free parameter. 

From the isomorphism (J6|, we have the K^ + ^(A)-matrices evaluated at 
A = given as, 
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-2A1) 
jii 
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4to,o)(0) = 
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> <J,o) (0) = 





An 















q 2 ) 




I 





q 2 J 



( 



K, 



(+) 
(0,0,1) 



(0) 



q- 2 
1 








(19) 



where 



A'l 



A3 



V 9 2 i5 y 



2 + (1 + g 2 )a sinh 7 
2 — g _2 asinh7 ' 
2 + (1 + q~ 2 )a sinh 7 
2 — g 2 asinh7 



A22 



2 + (g 2 + g 2 )asinh7 
2 — a sinh 7 



(20) 
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and 



41d (o) 



41o)(°) 



where 



/ q -2 








.[77] 

Jll 



\ 






,2 AH] 



o <rjii J y 
o \ 



iff] 



g 2 y 



a; 



(+) 

(1,0,1) 



(0) 








1 







o <fj& 



.[77] 

J33 



2 + q~ 2 a sinh7 
2 — (1 + g 2 )asinh7' 
2 + q 2 a sinh 7 



.[77] 
J22 



2 + a sinh 7 



2 — (q 2 + q 2 )asinh7 



2 — (1 + g _2 )asinh7' 
ight boundary free parameter. 



(22) 



and a is the rij 

Therefore, using (16) we can write the corresponding right boundary 
terms 



B 



( 2—q 2 asinh7 
4 cosh 7 






2—q 2 osinh7 
4 cosh 7 



B 



(Hi]) 



V 

( 2— a sinh 7 
4 cosh 7 








2+(l+ij 2 )asinh7 
4 cosh 7 





2+(q 2 +q~ 2 )a sinh7 



V 







4 cosh 7 





B 



(Hi]) 



/ 2+(l+g _2 )asinh7 
4 cosh 7 








2—q 2 a sinh 7 
4 cosh 7 








2— o sinh 7 
4 cosh 7 






2—q 2 a sinh7 
4 cosh 7 / 



(23) 
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and 



B 



(!.[«]) 



B 



(2,[/J]) 



5 



(3,[/i]) 



' 2+g 2 asinh7 
4 cosh 7 






V 

f 2+a sinh 7 
4 cosh 7 








2+q~ 2 a sinh7 
4 cosh 7 








2— (l+q 2 )a sinh 7 
4 cosh 7 



\ 





2—(q 2 +q~ 2 )a sinh7 
4 cosh 7 



V 

/ 2~(l+g~ 2 )asinh7 
4 cosh 7 













V 







2+g 2 asinh7 
4 cosh 7 








2+a sinh 7 . 
4 cosh 7 / 




2+q 2 a sinh 7 
4 cosh 7 



/ 



(24) 



are allowed B^'^'^' 



Here we have 6 different integrable boundaries related by the isomorphism 
6|). However, it is worth to note that other combination of the boundaries 

with i,j = 1,2,3, a,b — I, II re- 
sulting in 36 integrable boundaries for the spin-1 U q [sl(2)] Temperley-Lieb 
Hamiltonian. 

The action of the boundary terms on the Hilbert space is given by 



r>(i,j,[ a M) 



where 8^'^ 



Lacr T~ I TT 



^=S^ b]) \a---r> (25) 
and the sites are indexed by a, t = (1, 2, 3) = 



). Here we recall that laa and Ttt' 



are the matrix elements of the 



boundary matrices B±'^ and B K [ ,m> respectively. 

In the next section, we will restrict ourselves to the case of integrable 



boundaries related by the isomorphism (.B^^ = B\ l,[a>> + B { l' m> ) and we 
shall use the coordinate Bethe ansatz in order to obtain the eigenvalues of 
the Hamiltonian 081). 



»(i.W) 



>(i,[a])N 
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3 Coordinate Bethe ansatz 



In most of the cases where Bethe ansatz is successfully applied, one can build 
up all the eigenstates from just one reference state and usually there exist 
only few of such reference states available. In the case of Temperley-Lieb 
spin chains, by contrast, we have exponentially degenerated ground states, 
which implies that we have a very large number of reference states. In fact 
we have 3 x 2 L_1 natural eigenstates which can be used as reference states. 
This explains the difficulties in constructing all the eigenstates from just one 
reference state [9| [T0\ [TT]. 



£ u < 


1+ + ++) 
l+o + +) 

I+ + 0+) 


£12 < 


I+ + +0) 

l+o + o) 
l+o-o) 


£13 < 


V + o-) 
J+o--) 


£21 < 


✓ 

|0 + ++) 
|0 + 0+) 
|0 - 0+) 


£22 < 


|0 + +0) 

J0--0) 


£23 < 


f 

|o- ■— ) 

|0 — 0— > 
|0 + 0-) 

> 


£31 < 


'1-0++) 

J--o+) 


£32 < 


f 

1 o) 

|-0-0) 

l-o + o) 

1. 


£3 3 < 


1 ) 

I--0-) 
l-o--) 



Table 1: The reference states (natural eigenstates) of the Hamiltonian ([8]) 
and its eigenvalues for L — 4. 

The reason for such differences is that the bulk part of the Hamiltonian 
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p(i>[<*]) c(i>M) 


p(l,[a]) 
c 13 


p(i.N) _ _ _ c(i,N) 

°21 ~~ °22 ~~ °31 ~ °32 


_ £.(l,[o]) 
°33 ~~ °23 


p(2,[a]) p(2,[o]) 
^12 — ^32 


p(2,W) 
c 22 


c(2,H) t»(2,[a]) p (2,[o]) p (2,[a]) 
c ll — c 13 ~~ c 31 ~~ °33 


M2,[a\) ~(2,[a]) 
c 21 — c 23 


c(3,[a]) 

^11 — ^21 


p(3,H) 
c 31 


c( 3 >[°]) p( 3 >[°]) p(3,[«]) p(3,[o]) 
c 12 — °13 ~~ c 22 — °23 


c.(3,[a]) p (3,[a]) 
c 32 — c 33 



Table 2: The relation among the boundary eigenvalues for different solution 
of the reflection equation. This relations hold true for any a = I, II. 

(the Temperley-Lieb operator Uk,k+i) is the projector operator onto the two- 
site spin zero singlet. This implies that there exist 3 x 2 i_1 states which are 
eigenstates of the bulk Hamiltonian with zero eigenvalues. Therefore, these 
states are also eigenstates of the boundary part of the Hamiltonian B^'^ 
with eigenvalues Sar"^ (see e.g Table El). 

Moreover, apart from the natural degenerescence of the boundary eigen- 
values Ear , one can see from the structure of the boundary matrix 
that not all £ar are independent. In fact they are also degenerated and 
can be grouped in four blocks for each integrable boundary related by the 
isomorphism (see Table [2]). 

In face of the large number of reference states, the standard construction 
of the all eigenstates seems to be impracticable, though it is possible. How- 
ever, in order to obtain the eigenvalues of the Hamiltonian it is enough to 
work out with a few reference states. In fact, we can take one reference state 
from each block of eigenvalues . From now on, we drop the label for 
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different solutions of the reflection equation from the boundary eigenvalues, 

SUch that Ear"^ = Ear- 



3.1 Ferromagnetic reference state 

We shall start by considering the pseudo particle as a singlet over the stan- 
dard ferromagnetic state. Therefore, it is convenient to start our ansatz with 
the following linear combination of the basis states [TT] . 



|n(*)> = <f 



+ - + ••• + )-<? 1 



+ • • • + 00 + • • • + 



+ - + •••+,. 

(26) 



which is an eigenstate of U kk+ \ such that 



U kM1 \Q(k)) = Q \Q(k)) , U k+1M2 \Q(k)) = \Q(k + 1)) , 

U k>k+1 \Q(k ± 1)) = \Q(k)) , E7 fc _ 1)fc \Q(k)) = \Q(k - 1)) , (27) 

U k , k+1 \Q(j)) = 0, iik^{j,j + l}, 

where Q = (q 2 + 1 + q~ 2 )- Therefore, the action of the Hamiltonian H = 
Yl k Zi U k>k+ i + Bi L over this state results, 

H \n{k)) = Q \Q(k)) + £ u \tt(k)) + \n(k - 1)) + \Q{k + 1)) , 1< k < L - 1(28) 



h\q{i)) = Q|n(i)) + £u|n(i)) + |n(o)) + |n(2)) (29) 
h\q(l-i)) = Q\n(L-i)) + £i 1 \n(L-i)) + \n(L-2)) + \n(L)iqo) 

where B hL \Q(k)) = £ n |fi(Jfe)) ,l<k<L-l, |fi(0)> = {B 1>L - £ u ) 

and \Q(L)) = (Bn — £ n ) \Q(L — 1)). In addition to the previous relation, 
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we have a set of closing relations 

Him) = a! x) \m)+£v iA i (si) 

h\q(l)) = AW\n(L-i)} + s lu ^\n(L)}, (32) 

where A ; (1) = (£ 2l - £ n ) + q 2 (£ 31 - £ n ), A« = (£ l2 - £ u ) + q' 2 (£ 13 - £ n ). 
In order to cover all the solution of the reflection equation, we introduce the 
following notation v\ = (2, 2, 3) and u\ = (3, 2, 2) whose elements i — 1,2,3 
represent different solution of the reflection equation (for any a = I, II). In 



the above relations, we exploited the fact that, 

tfi >2 |n(0)> = A z (1) |^(l)) , (33) 

B hL \n(o)) = £ Vi:ll \n(o)), (34) 

U L ^ L \Q(L)) = AW\Q(L-1)), (35) 

B ltL \Sl(L)) = £ 1Uitl \n(L)). (36) 



3.1.1 One-particle state 

In the first non-trivial sector, we assume the following ansatz for the eigen- 

states 

\^) = J2Mk)\m)- (37) 
k=i 

Imposing the eigenvalue equation H l^) = E 1 \^fi) is fulfilled, we obtain 
a set of equations for the function A(k), 

(Q + £ u - E 1 )A{k) + A(k - 1) + A(k + 1) = 0, 1< k < L - 1 (38) 
(£ v . ll -E 1 )A(0) + A ( i 1) A(l) = 0, (39) 
(S^-E^A^ + ^AiL-l) = 0. (40) 
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Taking the ansatz for the amplitude 

A(k) = a(6)i k - a{-6)C k , (41) 



and substituting in Eq.(38) provides the following expression for the energy 
eigenvalues 

E 1 = Sn + Q + Z + C 1 - (42) 
The parameter £ and the ratio of the amplitudes a(8) / a(—6) are fixed by the 



equations (39 40), which results in the Bethe ansatz equation 

Q+s n - s Vil i + r 1 + e(i - fQ+£u- £x Uil + r 1 + ea - 



^2L 



Q+£n- s Vi xl + e + - A f 1} ) / w + ^ii - ^ i + £ + £- '(! - 

(43) 



3.1.2 Two-particle state 

In the next particle sector, we have two interacting pseudo-particles, which 
can be represented as a product of two pseudo-particles eigenstates, as given 
by 

|tf 2 >= Mki,k 2 )\n(h,k 2 )), (44) 

fci+Kfca 

where 

i 

\Q(k u k 2 )) = {-l) i+j q i+j - 2 \H-i,i)-M-3,3)) ■ (45) 

We can split the action of the Hamiltonian on the state \fl(k\, k 2 )) in four 
cases: 

(i) The case where two pseudo-particles are separated in the bulk, 

H \tt(h, fc 2 )> = (2Q + £ u ) \Q(k lt k 2 )) + \Q(ki - 1, k 2 )) + + 1, k 2 )) (46) 
+ |0(^i, A; 2 - 1)) + \n(h, k 2 + l)), 1 <h + 2<k 2 <L-l 
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(ii) The case where the pseudo-particles are separated but one of them or 
both are at the boundaries 

H fc 2 )> = (2Q + S u ) fc 2 )> + |ft(0, fc 2 )) + |fl(2, fc 2 )> 
+ |Q(1, A; 2 - 1)> + |Q(1, A; 2 + 1)> , 5 < k 2 < L - 1 (47) 
H \n(h, L — 1)) — (2Q + 5n) L — 1)) + Mh — 1,L — 1)> 
+ ^(jfei + 1, L - 1)) + |fi(fci,L - 1)) + L)> , 1< fci < L - 3 (48) 

if |Q(1, L - 1)) = (2Q + £ n ) L - 1)) + |ft(0, L — 1)) + |n(2, L - 1)) 
+ |fi(l,L-l)) + |Q(l,f)>, (49) 

(iii) The case where the particles are neighbours in the bulk 

H \Q(k, k + 2)) = (2Q + S n ) \Sl(k, k + 2)} + \Q(k -l,k + 2)) + \Q(k, k + 3)) 
+ \Q(k + l,k + 2)) + \Q(k,k + l)) , l<k<L-3 (50) 

(iv) The case where the particles are neighbours at the boundaries 

H\Q{1,3)) = (2Q + £ U ) 1^(1,3)) + 1^(0,3)) + |0(2, 3)> 
+ |0(1,2)) + |0(1,4)> , (51) 

H\Q(L-3,L-1)) = (2Q + £ 11 )\Q(L-3,L-l)) + \n(L-A,L-l)) 
+ \Q(L - 2, L - 1)) + \Q(L — 3, L — 2)) + \Q(L - 3, L)) . (52) 

In the above relations, we have introduced new states whose definition 
are given by 

\Q(0,k 2 )) = (B hL -£ 11 )\Q(l,k 2 )), (53) 
\n(h, L)) = (B hL - S n ) Mh, L — l)), (54) 
\n(0, L-l)) + |0(1, L)) = (B hL - S n ) Ml, L — l)) , (55) 

\n(k + i,k + 2)) + \n(k, k + i)) = u k+lM2 \n(k, k + 2)). (56) 
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The action of the Hamiltonian on these states can be written as follows 

H |fi(0, k 2 )) = A[ 1} fc 2 )) + £ Vi ll |fi(0, fc 2 )) + |fi(0, A; 2 - 1)) 

+ g|fi(0,A; 2 )} + |fi(0,fc2 + l)), (57) 

# in^i, l)) = aw l - 1)) + £ lMM £)) + - 1, l)) 

+ Q + \n(h + l,L)), (58) 

if |n(Ar, fc + 1)} = Q \Q(k, k+l)) + \Q(k - 1, k + 1)) + |fi(Jfe, fc + 2)) 

+ £ U |0(A;, & + !))■ (59) 

In order to obtain the eigenvalues, we have to substitute the above rela- 
tions in the eigenvalue equation (H |\I/ 2 ) = E 2 l^))- This will provide us the 
following set of equations for the amplitude A(k\, k 2 ), 

(2Q + S n - E 2 )A(k h k 2 ) + A(ki - 1, k 2 ) + A{k x + 1, k 2 ) + A(k h k 2 - 1) 
+ A(k h k 2 + 1) = 0, (60) 
(Q + S u - E 2 )A(k, k + 1) + A(k - 1, A; + 1) + A(k, k + 2) = 0, (61) 

(Q + - £ 2 )A(0, fc 2 ) + A(0, k 2 - 1) + A(0, A; 2 + 1) + A[ 1) A(1, fc 2 ) = 0, 

(62) 

(Q + ^ - E 2 )A(ki, L) + A(A;i - 1, L) + A{h + 1,L) + A^A{h, L - 1) = 0. 

(63) 



One can obtain the eigenvalues from the equation (60), 

E 2 = £ u + 2Q + & + er 1 + 6 + (64) 

provided that the following parametrization for the amplitudes is assumed 

A(k u k 2 ) =Y,e P a(9 l ,9 2 )fr$', (65) 
p 
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where the sum extends over all permutations and negations of momenta (9j), 
such that £j = e l9i , and Ep is the signature of permutations and negations. 
This structure already reflects the existence of the boundary. 



On the other hand, the equation (61 ) is the meeting condition for the two 



pseudo-particle states. Using the ansatz (65), we obtain the following phase 
shifts, 



a(9 2 ,9 1 ) 

a{e 2 ,-e x ) 
a{-e 2 ,e x ) 
a{-e 2 ,-e x ) 



3(92,-0!) 

s(-9 u 9 2 ) 

s(-02,0l) 

s(9 1 ,-9 2 ) 
si-92,-9, 



s{-9i,-9 2 



a{9 u 9 2 ), 

a(-eM, 

a(9 1 ,-9 2 ), 
a(-9 u -9 2 ), 



where 



s(9 1 ,9 2 



(66) 
(67) 
(68) 
(69) 

(70) 



At this point, we still have two remaining equation (62 63) which intro- 



duce the boundary effects. One can introduce the expressions (64 65) in the 



equation for the left boundary (62), which results 



where 



Fi( 



a(-9 1 ,9 2 ) = F l {9 1 )a{9 1 ,9 2 ), 
a(-9 2 ,9 1 ) = F l (9 2 )a(9 2l 9 1 ), 



o+f u -p 1 +^ 1 +e 1 - 1 (i-A! 1) ) 



(71) 
(72) 



(73) 
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Likewise for right boundary, one obtain the the following relations 

a{e 2> -e 1 ) = ^F r {9 1 )a{e 2 ,e 1 ), (74) 
a(6 1 ,-e 2 ) = e 2 L F r (e 2 )a(d 1 ,d 2 ), (75) 

where 

f (e 1 ) = ( g + gn-^i + ^ + ^d-A^ N 



Combining these relations with the phase shift relations (66 69 ), we obtain 
the Bethe ansatz equations 



^^)W>(|^)(f^j. (78) 
3.1.3 General sector 

The generalization to any number n of pseudo-particles goes along the same 
lines as before, although the calculation becomes cumbersome. Therefore, 
we just present the final results. The energy eigenvalues are given by the 
sum of single pseudo-particle energies 

n 

E n = £ u + J2Q + ^ + &\ (79) 
i=i 

and the corresponding Bethe ansatz equations depend on the phase shift of 
two pseudo-particles and on the boundary factors 

3=1 
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3.2 Other reference states 

In order to obtain the whole spectrum of the Hamiltonian we have to consider 
additional reference states. This has to be done for each different boundary 
eigenvalues £ CT . As a result of that, we must have as many as reference states 
and consequently Bethe ansatz equations as boundary eigenvalues. 

In principle, we have nine boundary eigenvalues £ ar . If one choose one 
reference state for each boundary eigenvalues (e.g the first state of each block 
of Table [T] extended to L-sites) and proceed along the same lines as the 
previous section, we obtain nine eigenvalues expressions 



it 




(81) 



i=l 



as well as its associated Bethe ansatz equations 




(82) 



where 




Q + SaT-SvurT+^ + Ml-^) 

Q+s aT - s aUiT + a+ er\i - ^ 
Q + s aT - £ aUiT + er 1 + 6(i - a^ t) ) 




(83) 
(84) 




(£31 - £2i)q 2 + (£11 - £2i)q~\ 
{£11 ~ £u)q 2 + {£13 ~ £i2)q~\ 



(85) 




(£21 - £31) + (£11 - £31) 9 



-2 



(86) 



(87) 



A (3) 

r 



(£ n -£i3)g 2 + (£i2-£i3). 



(88) 
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The remaining: index defined by v 2 — (1, 1; 3), v s = (1, 1> 2) and the 

u i>T are given by m 2 = (3, 1, 1), u 3 = (2, 1, 1). 

However, we can see from Table |2]that most of these equations degenerate 
into each other, resulting in four equations for each integrable boundary. We 
have verified numerically the completeness of the spectrum up to L = 6 sites. 

4 Conclusion 

In this paper we obtained the spectrum of the spin-1 U q [sl(2)\ Temperley- 
Lieb spin chain with diagonal open boundary conditions. We have identified 
that this model has large number of possible reference states. By selecting a 
small subset of these states, we manage to obtain four eigenvalue expressions 
and its associated Bethe ansatz equations by means of a generalization of 
the coordinate Bethe ansatz. This provides the complete description of the 
spectrum of the model for any values of the boundary parameters. We verified 
that the Bethe ansatz results are in agreement with the direct diagonalization 
of the Hamiltonian up to six sites. 

Apart from the news results, we believe that this work also brings a 
better understanding on the coordinate Bethe ansatz construction of the 
eigenstates. Although with this new perspective is possible to construct all 
the eigenstates for finite system size, this seems to be rather impracticable. 
Therefore, we still leave the problem of counting of the spectral multiplicities 
as an open question. We also hope that this work would shed some light on 
the algebraic Bethe ansatz construction for the Temperley-Lieb spin chains. 
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